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The recent PRB 73, 224501 (2006) henceforth referred as Ref.Q asserts that self-heating decreases 
with sample area reduction and claims to identify the intrinsic cause of ITS in submicrometre 'mesa'. 
I will show that this assertion lacks substantiation. I will further demonstrate that one and the same 
R{T) and the parameter-free Newton's Law of Cooling describe quantitatively a rich variety of ITS 
behaviours taken by Ref. [l| above and below T c at bath temperatures spanned over 150K. Thus this 
finding presents strong evidence in favour of heating as the cause of the 'intrinsic tunnelling spectra' 
(ITS) promoted by Ref.[J|. 



PACS numbers: 74.45.+C, 74.50.+r, 74.72.-h, 74.25.Fy, 

Heating is probably the most common problem in low 
temperature research and a harsh limiting factor for the 
study of current-voltage characteristics (IVC). Its signa- 
tures are particularly well studied for the case of met- 
als and conventional superconductors; notably, heating 
is known to cause IVC nonlincarities and transform a 
single- valued IVC into a multi-valued characteristic with 
regularly spaced branches (see Ref. [2| for review). Un- 
like conventional superconductors, high temperature su- 
perconductors (HTSC) reveal exceptionally poor thermal 
and electrical conductivities which makes them particu- 
larly prone to local heating. The heating-induced IVC 
nonlincarities in such circumstances were found to re- 
semble true tunnelling characteristics so closely that one 
may easily fall into the trap of identifying HTSC mate- 
rial as an SIS Josephson junction (see Ref.[3[ for review). 
The IVC by the authors of |lj reproduced in Fig. 1(a) 
are remarkably similar to extrinsic ones; however they 
are argued to be virtually free of heating. Below I will 
address the consistency of the basic assumptions, the ex- 
perimental data and the conclusions by the authors of 

h 

Assuming that self-heating in samples of different area 
is proportional to heat W=IV only, the authors of Ref. [l| 
claim that self-heating can be significantly reduced by 
means of sample area reduction. However, this assump- 
tion is incorrect as heat W, dissipated in a sample, es- 
capes through its surface area A, so the temperature rise 
depends on the heat load P — W/A. Hence this claim 
lacks grounds and, moreover, is at odds with experimen- 
tal data which point to the area independence in heating 
effects, see Ref.[3| for review. Additional evidence in sup- 
port of this conclusion may be seen in Ref. 4] , which finds 
that practically the same heat loads P ~ 10kW/cm 2 
build the ITS gap in the mesas of vastly different area 
1 < A < 30fim 2 made of the same Bi2212 crystal. [ll| 

As far as temperature rise is concerned .the systematic 
experimental studies summarised in Ref.Q suggest that 
the mean temperature, T, of the self heated sample is 



appropriately described by Newton's Law of Cooling, 



T = T B 



P/h, 



(1) 



where Tb is the temperature of the coolant medium (liq- 
uid or gas) and h is the heat transfer coefficient, which 
depends neither on A nor T . Furthermore, [3[ shows that 
in layered HTSC heating-induced IVC nonlinearities ex- 
ceed the intrinsic ones so radically that the latter might 
be safely ignored even at quite modest overheatings. This 
finding, also supported by the experimental data in [l[, 
is of particular importance to the present discussion. 

Although the assumptions by the authors of Ref. [l[ are 
not beyond dispute, it is worth considering their data, 
obtained in a state-of-art experiment. First I will verify 
the experimental consistency of [l[ using the IVC ori- 
gins least affected by heating and thus allow comparison 
with R(T) data measured independently. As shown else- 
where, Q, in the absence of heating the initial slope of 
IVC taken at Tb corresponds to normal state resistance, 
R N (T B ), which merges R(T) if T B > T c . Indeed, close 
correlation between the thus determined Rn(Tb) and 
R(T) measured by [l[ above T c is clear from Fig. 1(b). 
This correlation provides strong evidence in support of 
the experimental consistency of the data by [l[ . Further- 
more, in qualitative agreement with the direct measure- 
ments by H, 0], Rn(Tb) in Fig. 1(b) continues its growth 
and retains its upward curvature when the temperature 
is lowered through and below T c , as shown by the thin 
line in Fig. 1(b). However, as seen in Fig. 1(b) deviations 
from this behaviour rapidly develop when Tb lowers be- 
low 60K. As shown by Yasuda et al,0] heating is a likely 
cause of such deviations, so there is a possibility that the 
lowest two Tb were seriously underestimated (see the ar- 
row in Fig. 1(b)). Thus our analysis demonstrates the 
reasonable consistency of the experimental data by [l[ 
and suggests that in at least five of seven IVC by [l| the 
quoted Tb might be used as a reliable starting tempera- 
ture in the analysis of heating issues. 
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FIG. 1: (a): Solid lines reproduce nonlinear IVC, reported by 
Ref.[l[ for different Tb=4. 2,55, 75,85,95, 130.230K above and 
below T c — S9K; broken lines show the levels of constant heat 
load P=IV/A=5,10,20,50,100kW/cm 2 ; A^QM^m 2 . (b) 
compares measured R(T) shown by the thick solid line and 
R(P— >0) vs Tb estimated from the initial slopes of the IVC 
in Fig. la. 



Thus, the experimental data by the authors of [l| are 
sufficiently consistent and make it possible to address the 
origin of IVC nonlinearities using the parameter-free de- 
scription by Rcf.[3]. These IVCs, which are central to 
the discussion in are reproduced in Fig. 1(a) together 
with the levels of constant heat load. As could be eas- 
ily seen from this figure, in remarkable similarity with 
the basic ITS studies, nearly lOkW/cm 2 is required to 
build the characteristic IVC features attributed to the 
superconducting gap (ITS gap) by [lj. Such loads ex- 
ceed the critical ones by several orders, hence suggesting 
strong heating. However, additional analysis is required 
to quantify the temperature rise and, most importantly, 
to verify the extent to which the IVC nonlinearities might 
reflect the plausible intrinsic ones. 

To discriminate between intrinsic and extrinsic contri- 
butions, let us compare the measured data with those 
calculated on the basis of the assumption that heating- 
induced IVC nonlinearities exceed the intrinsic ones so 
radically that the latter might be safely ignored. The 
IVC in such circumstances is primarily determined by 
i?Ar(T), while the self heating is appropriately described 
by Eq.(l), see [3| for details. To make such quantitative 
analysis feasible, the heat transfer coefficient, h should 
be determined for the experimental conditions of [l[ . As 
shown in |3| (and reaffirmed by independent measure- 
ments by [a]), h could be obtained from a single self- 
heated IVC provided that the sample's R(T) is known. 
As will be shown in the next paragraph, the good data 
presentation by the authors of [l[ makes it possible to 
fulfil this task with even higher reliability. 



FIG. 2: (a): Compares measured R(T) shown by the thick 
solid line with those calculated with Eq.(l) from five of seven 
nonlinear IVCs from Fig. 1(a) using the same heat trans- 
fer coefficient h — 350W cm~ 2 K~ 1 for the data taken at 
Tb spanned over 150K; solid dots represent corresponding 
R(P->0) vs T B = 75, 85, 95, 130, 230K. (b): Solid lines re- 
produce the IVC from Fig. 1(a); broken lines show the lev- 
els of constant overheating AT=T - T s =15,30,60,150,300-Jf; 
A=0.36fim 2 . 



In terms of heating issues, it is appropriate to consider 
R = V/I as a function of heat load, P = VI /A, rather 
than IVC only (see above and also Ref.0). Provided 
that Tb is correct, the thereby determined R(P) could 
be converted with Eq.(l) into R(T) thus giving a reli- 
able estimate of h. The set of drastically different IVC of 
the same sample at various Tb reported by [l| provides a 
harsh consistency check for our approach as one and the 
same h should convert various IVC into a single R(T) 
thus leaving absolutely no space for manoeuvre. How- 
ever, as seen from Fig. 2(a), the parameter-free Eq.(l) 
collapses all IVCs obtained at Tb spanned over 150K into 
a single curve which reproduces quantitatively the R(T) 
of the same 'mesa' and allows estimate of the heat trans- 
fer coefficient /i= (300-350) Wcm" 2 ?^ 1 . Thus, Fig.2(a) 
confirms the heating origin of the IVC non-linearity and 
suggests that the IVC by Ref.[l[ will be almost linear 
above and below T c = 89K if the heating artefacts are 
removed. 

In addition it is worth considering the remaining two 
IVCs measured in Ref.[l| in this sample. Albeit there 
are no physical reasons to believe that these data are 
unaffected by heating, it appears that for Tb=4.2K, the 
authors claim that 'no values of h can be found to fit the 
experimental R(T) curve satisfactorily even if the specific 
shoulder structure in the IVC were not present'. The key 
to the resolution of this confusion is presented by Fig. 1(b) 
and the corresponding discussion which suggest that the 
quoted values are seriously underestimated as compared 
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to the effective Tb (see the arrow in Fig. 1(b)). The extent 
of this underestimation could now be evaluated, since the 
very same h = 350W ~cm~ 2 K~ x should convert these IVC 
into the very same R(T) hence providing an independent 
way to estimate the effective Tb- Our analysis suggests 
that these IVCs are actually taken at Tb —60 and 65K 
correspondingly, thus reaffirming quantitatively the val- 
ues anticipated in Fig. 1(b). 

Finally, using this h—SBOWcm^K -1 one can easily 
quantify the self-heating which builds any IVC point of 
interest (see Fig. 2(b)). The data in this figure suggest 
that the temperature rise calculated in Ref.[lJ under the 
dubious assumption that the heat dissipated in the sam- 
ple, sandwiched between metal electrode and substrate 
of exceptionally poor thermal conductivity, escapes ex- 
clusively into the substrate is seriously underestimated. 

To conclude, it is demonstrated, using exclusively the 
data from the commented article, that unlike the remark- 
ably consistent state-of-art experiment by [l| , neither the 
interpretation nor the conclusions are beyond dispute. It 
is shown that the experimental IVC taken above and be- 
low T c at vastly different Tb spanned over 150K are de- 
scribed quantitatively by Newton's Law of Cooling and 
Ohm's law using the normal state resistance of the same 
sample only. This finding confirms the heating origin 



of the IVC by_ll| and suggests that unlike conventional 
spectroscopy [9(, the heating in ITS is not a small per- 
turbation but a principal cause of IVC nonlinearity, no 
matter whether the sample is of centimetre or submi- 
crometre size. 

Our conclusions by no means rule out the experi- 
mental approach by the authors of [l|. In addition to 
the remarkable consistency mentioned above, this ap- 
proach provides a greatly improved heat transfer coef- 
ficient which makes it possible reliably to address the 
intrinsic response and perform other worthwhile IJT ex- 
periments, some of which were proposed by Indeed, 
the h = 300 — 350Wcm~ 2 K _1 estimated above repre- 
sents a very major improvement, as so far all known ex- 
perimental h fell into the 2-60Wcm~ 2 K~ 1 range 
Furthermore, this improvement seems to be robust as 
the quantitatively similar h=30QW 'cmT 2 K^ 1 is revealed 
by another sample (of somewhat smaller A = 0.09/im 2 ) 
by this group, [10l |. 
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